Two-dimensional relativistic hydrogenic atoms: A complete set of constants of motion 
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The complete set of operators commuting with the Dirac Hamiltonian and exact analytic solution 
of the Dirac equation for the two-dimensional Coulomb potential is presented. Beyond the eigenvalue 
^ of the operator jz, two quantum numbers ri and k are introduced as eigenvalues of hermitian 
operators P = /3a'^ and K — f3{a'^lz + 1/2), respectively. The classification of states according to the 
full set of constants of motion without referring to the non-relativistic limit is proposed. The linear 
Paschen-Back effect is analyzed using exact field-free wave-functions as a zero-order approximation. 
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I. INTRODUCTION 

Low-dimensional quantum systems have been the focus 
of extensive theoretical investigations in the last decades. 
Technological advances in semiconductor physics and re- 
cent developments in nanostructure technology provide 
techniques of creating low-dimensional structures like su- 
perlatices, quantum dots, quantum wires or quantum 
wells jH]. The most representative analogues of hydro- 
genlike systems in the world of semiconductors are hydro- 
genic donors being the bound states of conduction elec- 
tron and a donor impurity Q and Wannier-Mott excitons 
formed by an electron and a hole . After the renormal- 
ization of Coulomb potential by introducing of dielectric 
constant and replacing the electron mass by effective one, 
the atomic objects in two-dimensional structures can be 
treated as 2-D hydrogelike atoms. At this stage, the 2- 
D hydrogen problem determines a leading approxima- 
tion for study of hydrogen type bound states in extreme 
anisotropic crystals, in which the z-component of a di- 
agonal anisotropic mass tensor is much larger then the 
two remaining ones [4|. The non- and weak-relativistic 
approaches are usually considered as sufficiently good 
approximations to the realistic description of the 2-D ob- 
jects in solid matter. However, in a searching of the 
quantum-mechanical properties of the 2-D systems in- 
teracting with low-dimensional gauge fields [1, 0|, the 
complete relativistic theory is inevitable. 

The quantum mechanical two-dimensional central 
problem, with the Coulomb potential —Z/p has been 
solved by many authors. In the nonrelativistic theory 
the analytic solution can be derived in strict analogy to 
the three-dimensional Coulomb problem, after the sepa- 
ration of the Schrodinger equation in polar coordinates 
or in parabolic coordinates [lOj. The solution of the 
relativistic 2-D hydrogen-like problem has been obtained 
in the framework of the two-component approach fll[. 
Although the formalism based on the two-dimensional 
representation of gamma matrices gives correct formula 
for energy levels, it does not provide a good background 



for analysis other observables. Alternatively, standard 
Dirac-Pauli representation of the Dirac matrices can be 
used. The four-component analytic solution of the Dirac 
equation with the Hamiltonian 



H = ca-Tp + /3c 



Z 
P 



(1) 



in two spatial dimensions has been obtained by Guo et 
al. [T^ . The authors have investigated two decoupled 
eigenstates of H 



/i(p)e'(^-i/2)<^ 



i5i(p)e»(^+i/2)0 





2.g2(p)e'(f-i/2)0 




(2) 
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where /i is the eigenvalue of the z— component of the total 
angular momentum jz defined as: 



I. 



with 



and 



(Jz 

a. 



-idld(f). 



(4) 



(5) 
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In non-relativistic limit the states ^ and ([3]) describe an 
electron with spin wp and down, respectively. Without 
any additional information about other conserved quan- 
tities the radial amplitudes of these two decoupled states 
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have to be determined by two different sets of radial equa- 
tions [l^l- It seems a little undesirable, in particular if 
we appeal to the three-dimensional theory where only one 
pair of radial equations appears for all quantum states. 

The purpose of this paper is the analysis of integrals of 
motion of the two-dimensional relativistic hydrogen atom 
and classification of the states based on good quantum 
numbers. These goals are achieved by introducing, into 
Dirac Hamiltonian, hermitian operators associated with 
conserved quantities. Another important aspect of the 
presented approach is the possibility of determining all 
radial functions from one system of radial equations. 



II. EXACT SOLUTION FOR FIELD-FREE 
ATOM 

We have find, that two linearly-independent states ([2|) 
and ([3]) are the eigenstates of an operator 



(7) 



which commute with the Hamiltonian ^ . It follows from 
([7]) that P is an involution {P^ — 1) and has two eigen- 
values ?7 = ±1. In the non- relativistic limit the different 
signs of r/ correspond to states with opposite spin direc- 
tions. Moreover, it can be directly verified, by computing 
relevant commutators, that beyond operators jz and P 
there exists an operator K 



K = p{a'Jz + -), 



(8) 



which commutes with H and both operators P, jz ■ The 
eigenvalue of K can be referred to the Dirac quantum 
number k. To determine a physical meaning of quantum 
number k, let us derive a relation between k and [i. If 
we consider the square of K, 



K^^{a%f+a% + \ = {lz + \a'zf, (9) 



we obtain 



(10) 



which means that eigenvalues of K satisfy the relation 

K = ±\^i\. (11) 

We note that operators K, P and jz are not independent. 
They fulfil the relation 

K - Pjz. (12) 
which gives a similar relation for eigenvalues 

K = jirj. (13) 



It appears from ()13p that the upper and lower signs in 
(llip distinguish between two different physical situations. 



when, in non-relativistic limit, spin is parallel or antipar- 
allel to the total angular momentum. 

Now we are in a position to solve the 2-D Coulomb 
problem in strict analogy to 3-D case and perform the 
classification of states free of non-relativistic quantum 
numbers. To this end we introduce to the Dirac equation 
quantum numbers associated with the complete set of 
commuting operators {H^K,jz). 

In atomic units and polar coordinates {p, 0) Hamilto- 
nian ([1]) can be written in the form 



(14) 





H = c{apPp ^ 


- -apa'zlz) - 


P 


where 










(T-p 
(T-p 


, cr p = 


e-"f' ' 
e"*' 


and 




. d 





Pp 



dp' 



(15) 



(16) 



Taking into account the algebraic properties of matrices 
a, (3 and definition ([8|) of if we obtain 



H — dapTTp + -ap(3K) + [3c^ , 

P P 

where the radial momentum tt^ is defined as 

.,d 1 , 



(17) 



(18) 



In the representation in which operators H, jz and K 
are diagonal, energy levels are determined by radial part 
of the wave function only. The pertinent radial Dirac 
equation takes the form 



[ciapTTp + —apl3) + (3c^ ]R = WR, (19) 



where W denotes energy and Z is the charge of the 
Coulomb field. Since ttp and /3 fulfil relations 

a2:=/32 = l, apl3 + l3ap = Q (20) 

they can be represented by two dimensional hermitian 
matrices 



-i 

1 



1 

-1 



(21) 



According to Eq. (PT|) the radial function R has two 
component, which for convenience we take in the form 



R{p) = 



1 



P 



1/2 



F{P) 
G{p)Zlc 



Introducing two new variables 



Zp, E - 



W 



^2 



(22) 



(23) 
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and substituting (|2T|) into (|19p leads to the wave equation 
for the electron moving in the two-dimensional Coulomb 
field in the form of the pair of K-dependent radial equa- 
tions: 



ar r r 



(24) 



dF K 1 

F-[A(-+S) + 2]G = 0, (25) 

ar r r 

where A = {Z/c)'^. 

Since asymptotic solutions of Eqs. ([24| . ((25|) decay ex- 
ponentially we try to find the solution of radial equations 
(l24t and (1251) in the form 



i=0 



^ ^ j,7g-"'-^a.r', G = r'^e""'' ^ 6,r\ (26) 
where 



^J~E{2 + \E). 



(27) 



Substituting expansions (^5)1 into equations (^5)1 we 
obtain the linear relations between the expansion coeffi- 
cients 



(i -|- 7 — K)ai — Xbi = aui-i + {2 + \E)bi 



(28) 



(29) 



Oi + (i + 7 + K)hi — —Eai-i + Q;6i-_i . 

For 1 = 0, we have 

(7 — K)ao - Xbo = 0, ao + (7 + i^)bo = 0. (30) 

Since both ag and 69 are different from zero, the secular 
determinant must vanish, which leads to 



7 



±\/k2-A, 



and 



60 



7 -I- k' 



(31) 



(32) 



For j > 0, the expanding coefficients can be calculated 
iteratively from the relations 



where 



(i + 7 — K + a/E)wi 
i{i + 27) 



Oi = - (« + 7 + 



Wi = -Ea,_i + abi^i 



(33) 



(34) 



(35) 



The condition of square integrability of the wave- function 
allows only the upper sign in Eq. pi|) and requires the 
termination of power series ()26|) at a some power n' 



which leads to the condition 

EUn' = abn' ■ 



(37) 



Dividing ((34| by ((33| for i — n' and taking into account 
(P7)) and ([57)1 leads to the equation for E 



^/~E{2 + \E){n' + j) = l + \E. 
Solving this equation and substituting 
7 = - A, 

we obtain 

-1/2-1]. 



^ A ^ (n' + y/^^2—x)2 ' 



(38) 



(39) 



(40) 



In order to compare the expression (|40p with non- 
relativistic one, we define the principal quantum number 
as follows 



n = n'+\ k \ +1/2. 



(41) 



Since > we must have \ k \< {n — 1/2). However for 
n' — 0, the number k must have positive value only. The 
absence of the k < for n' — follows from equations 
and p7p . which both imply 



(7 + k) > 0. 



(42) 



According to (|39p 7 is a real number smaller than | k \ 
and inequality (|42p can be satisfied only if k is positive. 
Therefore k must fulfils relation 



1 , 

\ K \<n-l. 

I 2 I- 

Finally, energy levels are given through 



E„ 



[(1 



A 



(n- I At I -1/2+\/k2-A)2 

In the non-relativistic limit (A — )■ 0) we obtain 

2 



(43) 

)-V2_i]. 

(44) 



E.„ 



(2n- 1)2 



(45) 



The complete spin-space-dependence of wave- functions 
is given by 



fnK{r)n^f,{(j)) 



where 



n = 1,2,3, 



(46) 



(47) 



K = 1/2, -1/2, 3/2, -3/2, (n- 1/2), (48) 



and 



a„'+i = 0, bn'+i = 0, a„/ 7^ 0, 6„' 7^ 0, (36) 



(49) 
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The cylindrical spinor rj„^((/i) is defined as 



k+A£„4(m-1/2)0 

^K+i£gi(Ai+l/2)0 
2/^ 



(50) 



and radial amplitudes have the form 

n-|K|-l/2 
1=0 



n-|K|-l/2 



— rl-i/'^p-ar 



E 



(51) 



(52) 



with coefficients determined by relations 

Alternatively, in a similar way as in 3-D case, the solu- 
tion of radial equations fM)) and P5|) may be expressed 
in terms of confluent hypergeometric functions 



F{r) = r''e-°n(K + -)Fi{r) - n'F2{r)], 
a 



G{r) 
where 



-rTe-"'-[(K + -)Fi{r) + n'F2(r)], 
a a 



i^i(r) = iFi(-n',27+l;2ar), 



-F2(r)-ii^i(l-n',27 + l;2ar). 



(53) 



(54) 



(55) 



(56) 



In order to introduce the classification scheme based 
on spectroscopic notation it is useful to define an orbital 
quantum number as 



l^K-- 

' 2 



(57) 



Note that this quantity according to satisfies the 
inequality / < [n — 1). 

In Table U we display, as an example, the lowest few 
states with principal quantum numbers n = 1,2,3. We 
can see that for a given n, the states with the same | 
K I are degenerate and the lowest energy corresponds to 
minimal value of | k |. According to for each value 
of K there are two possible values of /i. Therefore, the 
degree of degeneracy of the nn-th energy level is 2 for 
K = (n — 1/2) and 4 for | k |< (n — 1/2), respectively. 
It is worth to point out that, in a contrast to the three- 
dimensional case, due to the equality | k, \ — \ /i |, the 
states with | /i |<| k | do not occur. 



III. LINEAR PASCHEN-BACK EFFECT 



Coulomb field with charge Z and in a static uniform mag- 
netic field. In atomic units the relativistic Hamiltonian 
can be written in the form 



= ca • (p -t- A) -f (3c^ 



Z 
P' 



(58) 



Taking into account the standard four-dimensional Dirac- 
Pauli representation of the Dirac matrices and the vector 
potential A = B x p/2, for B = Bz perpendicular to the 
plane of transversal motion of the electron, we can write 



a - (p + A) 



ipPp 



Bp 

2 ' 



(59) 



Introducing into Hamiltonian H both operators P and 
K, defined in a previous section, we obtain 



H = c{apTTp + -UpPK + -BpapPP) + I3c^ 



Z 

p ■ 2 ^ p 

and appropriate radial Dirac equation in the form 



(60) 



in irjB 

[c{apTTp-\ api:i+—— 

p z 



panf3)+Pc^^-]R = WR. (61) 
P 



Substituting (PT|) into (pTjl and changing variables ac- 
cording to and and introducing a new variable 
B B/Z'^ leads to the wave equation, for the elec- 
tron moving in the superposition of the two-dimensional 
Coulomb field and constant homogeneous magnetic field, 
in the form of the pair of k- and 77-dependent radial equa- 
tions 



dG , K 
dr r 



\^Br)G + (i 
2 r 



E)F = 0, 



^-{:i + lr^Br)F-[X{-+E) + 2]G- 
dr r 2 r 



0. 



(62) 



(63) 



The complete spin-space description of eigenfunctions is 
the same as for the field-free atom. The only difference is 
in dependence of radial functions on both quantum num- 
bers K and fi. In consequence, magnetic energy shift may 
depend on the symmetry of the states. We investigate 
this problem in the first-order perturbation approach. In 
order to applying the perturbation formalism we rewrite 
Eqs. dMl) and ^ in 2 x 2 matrix form 



where 



(/i(°) B/i^i) - ES)<i> = 0, 



-1/2 ~d/dr-K/r 
d/dr-n/r -(2 + A/r) 





-r/r/2 



-r]r /2 




(64) 



(65) 



(66) 



Let us consider now the Dirac Hamiltonian describing 
transversal motion of an electron around a fixed center of 



1 
A 



(67) 
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$(r) 



Fir) 
G{r) 



(68) 



Perturbation expansions for energy and wave-function 

oo oo 

E = ^E^^B\ $ = ^$Wb* (69) 
lead to the following perturbation equations 

n 

/j(o)$(") _^ _ ^ £:W5'$("-*) ^ 0. (70) 

The zero-order equation (n = 0) 

(/i(o)-S(°)5)$(") =0 (71) 

is equivalent to the system of radial equations ([M)) and 
(P5)l for field-free atom. Under the condition of orthogo- 
nality 



($(0) c$W\ = 0, 



(72) 



for i > 0, the n-th order energy correction can be written 
in the form 



($(0)^/j(l)$(«- 



The calculation of the first-order energy correction 



{F, F) + A(G, G) 



(73) 



(74) 



can be performed in closed form using radial functions 
([SH]) and dm). For n' = (k =| /i |), we obtain 



i?« = £(27+l). 



(75) 



Taking into account that hypergeometric functions (j55p 
and l[56)) depend only on | k |, the dependence of func- 
tions F and G on k is due to the factor {k + 1/2) in front 
of Fi in ([53|) and jM])- It means that for n' > the 
first-order energy corrections have the general form 



-V- 



-A. 



where 



ad — he 
K^c^ — d? ' 



A. 



bd 



(76) 



(77) 



Radial integrals a,b,c,d depend only on and n' and 
are given in Appendix A. It appears from Eqs. (|75p and 
([77)1 that the k,- and /i-degeneracy of field-free levels is 
completely removed by external magnetic field. Table HIl 
lists the first-order energy shifts of levels corresponding 
to states with n = 1,2,3, presented in Table H] We can 
see that linear corrections E^^"^ essentially depend on the 
sign of both quantum numbers k and /i. Calculating the 
non-relativistic limit of ([75)1 and ([7S)l we obtain the values 



of energy corrections i?^'' = limA^o ^'•^^ which consist 
with the non-relativistic result 



= i(m + 2m,). 



(78) 



where m and nig are the eigenvalues of and — ^cr^, 
respectively. The non-relativistic formula (|78p may be 
derived from the Schrodinger-Pauli equation in a similar 
way as in the 3-D case (see for example (T^V 



IV. CONCLUDING REMARKS 

The exact solution of 2-D hydrogen problem which 
have been presented in this paper is consistent with fun- 
damental principles of quantum mechanics. Due to in- 
troducing good quantum numbers into wave equation, 
the problem of energy spectrum is solved exactly in close 
analogy to the 3-D case. Moreover, the wave-functions 
are classified according to the complete set of constants 
of motion and determined by one K-dependent system 
of radial equations. Although, the energy spectrum of 
field- free central problem depends only on | k | = | /i |, 
the full K- and /i-dependence appears when the external 
magnetic field is applied. 



APPENDIX A 

Radial integrals occurring in the expression of the first- 
order magnetic energy correction ([75)1 are defined as fol- 
lows 

a=^Xi, h^-[{K' + \)K^-n'^K2\, (Al) 
a a a 



2, A£;\ ,,\E'^ 
c = -(1 + -^)h + 'in'{— - l)/i2. 



(A2) 



d= (l + ^)[(«;' + A)/i+n''/2] 



2n'XE ^. 
+ — (^ - l)/i 



where 



/ r^'T+^e-^'^'-Fidr, 



I,. = 



where i = 1 , 2 and 



Ii2 - / r^'^e-^^'^Fif^adr. 
Jo 



(A3) 

(A4) 
(A5) 

(A6) 
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TABLE I: Relativistic quantum numbers, spectroscopic notation and energies for bound states with n — 1,2,3 for the two- 
dimensional hydrogen atom. The values of energy are computed with c — 137.03599976 13]. 



n 


n' =n- \ K.\ -1/2 


K 


l=\^-j 1 


Notation 


Energy 


1 





1/2 





lSl/2 


-2.000106514052 


2 


1 


1/2 





2Si/2 


-0.222234057055 


2 


1 


-1/2 


1 


2pi/2 


-0.222234057055 


2 





3/2 


1 


2p3/2 


-0.222223537086 


3 


2 


1/2 





3Si/2 


-0.080002897124 


3 


2 


-1/2 


1 


3pi/2 


-0.080002897124 


3 


1 


3/2 


1 


3P3/2 


-0.080000624824 


3 


1 


-3/2 


2 


3(^3/2 


-0.080000624824 


3 





5/2 


2 


3(^5/2 


-0.080000170405 



TABLE IL First-order magnetic corrections ij''^-' to the en- 
ergies for states with principal quantum numbers n = 1, 2, 3 
of the two-dimensional relativistic hydrogen atom. The upper 
and lower signs in the front of S^^^ are referred to the two val- 
ues of = ± I K I, respectively. The numbers in brackets are 
the powers of 10 by which the entries are multiplied. In the 
column 5 nonrelativistic values of linear magnetic corrections 
are given. 



n' 




State 









1/2 


lSl/2 


±0.49997337 


±0.5 




3/2 


2p3/2 


±0.99999112 


±1 




5/2 


3^5/2 


±1.49999467 


±1.5 


1 


1/2 


2Sl/2 


±0.49999704 


±0.5 




-1/2 


2Pl/2 


=F2.9586[-6] 







3/2 


3p3/2 


±0.99999680 


±1 




-3/2 


3^3/2 


±0.49999680 


±0.5 


2 


1/2 


3Sl/2 


±0.49999899 


±0.5 




-1/2 


3pi/2 


=Fl.0651[-6] 
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